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Abstract We present the first in-depth statistical survey of flare source heights 
observed by RHESSI. Flares were found using a flare-finding algorithm designed 
to search the 6-10 keV count-rate w hen RHESSI's full s ensitivity was available 



in order to find the smallest events ( Christe et all |2008). Between March 2002 



and March 2007, a total of 25,006 events were found. Source locations were 
determined in the 4-10 keV, 10-15 keV, and 15-30 keV energy ranges for each 
event. In order to extract the height distribution from the observed projected 
source positions, a forward-fit model was developed with an assumed source 
height distribution where height is measured from the photosphere. We find 
that the best flare height distribution is given by g(h) oc exp(— h/X) where A = 
6.1 ±0.3 Mm is the scale height. A power-law height distribution with a negative 
power-law index, 7 = 3.1 ± 0.1 is also consistent with the data. Interpreted as 
thermal loop top sources, these heights are compared to loops generated by a 
potential field model ( PFSS). The measured flare heights distribution are found 
to be much steeper than the potential field loop height distribution which may be 
a signature of the flare energization process. 



1. Introduction 

Since the launch of STEREO, there has been renewed interest in understand- 
ing the three dimensional structure of solar flares. Stereoscopic vision by the 
STEREO spacecraft, allow true source heights to be determined through paral- 
lax. Determining the height of solar flare sources is important in order to identify 
the location of the acceleration region and constrain flare models. In the standard 
flare model, the electron acceleration region is situated in the corona above and 
produces closed magnetic loops which fill with hot thermal plasma observable 
in hard X-rays (HXR) up to 30 keV or more in large flares. At the footpoints 
of these magnetic loops, accelerated particles are stopped by collisions and are 
observed down to 10 keV for the smallest flares. 
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Past studies have used a variety of different instruments and methods to 
determine the height of flare sources on the Sun. Using simultaneous observations 
of 607 flares by Explorer 33, 35, and Mariner V, ICatalano and van Alien! f|l973h 
derived the height of soft X-ray sources (SXR, 1-6 keV) at the limb through the 
difference in occultation between observations with different view angles. They 
found that their measurements of the dependence of emission on altitude above 
the photosphere was well fit by an exponential (e~' l / A ) with a scale height of 
A = 11 ± 3 Mm. iKane et all (| 19791 ) determined that a HXR source (>50 keV) 
was "well below" a height of —25 Mm using simultaneous observations by the 
International Sun Earth E xplorer 3 (ISEE-3) and the Pioneer Venus Orbiter 
(PVO). Later, Kanej (1983) showed that for impulsive HXR sources (150 keV) 
95% of e mission originated a t a height of <2.5 Mm above the photosphere. Using 



Hinotori, lOhki et al 



( 19831) studied two impulsive HXR (17-40 keV) bursts 
near the limb and fo und the height of the source centroid was <7 Mm. Later, 
Tsuneta et ail (|l984l) compared Hinotori and Ha images and found the height 
of a HXR (16-38 keV) source to be -40 Mm. iTakakura et all (|l986l) compared 
HXR and Ha positions and fo und a height of 7 ± 3.5 M m for HXR (20-40 keV) 
sources. Using Yohkoh/HXT, iMatsushita et all (jl992[ ) performed a statistical 
study comparing source locations in Ha and HXRs (14-23 keV). Interpreting 
the difference geometrically, she found an average source height of 9.7 ± 2.0 Mm 
and that the source height decreases as a function of energy. The difference 
in heights were found to be -1.0 ± 0.3 Mm, -2.0 ± 0.5 Mm, -3.2 ± 0.7 Mm 
in the four Yohkoh/HXT energy bands, 23-33 keV, 33-53 keV, 53-93 keV, re- 
specti vely, relative to the lower energy channel. lAschwanden. Schwartz, and Alt 
(|l995l) using CGRO/BATSE, investigated milliseconds time differences between 
peaks in emission in the 25-50 keV a nd the 50-100 keV energy band. Interpr eted 
as a difference in the time of flight, lAschwanden. Schwartz, and Altl (|1995|) de- 
rived an altitude of 7.3 ± 0.8 Mm for the acceleration site. Using RHESSI, 
Aschwanden. Brown, and Kont ar ( 2002) determined the absolute height of foot- 



points at different energies (from 15 to 50 keV) which occurred near the limb by 
assuming a thick-target model and a power-law density model. Source heights 
between 4. 5 Mm (at 15 keV) to 1 Mm (at 50 keV) were found. More recently, 
Sato et all (120061) re-analyzed the Yohkoh flare catalog using the method pio- 
a|Mj 



neered by IMatsushita et all (|1992l) . Comparing sources in the four Yohkoh/HXT 
energy bands, it was found that the difference in source heights between each 
energy band was 460 ± 40 k m, 880 ± 90 km, and 93 ± 150 km, respectively, 
smaller than those found by Matsushita et all dl992l). Recen tly, in two papers 
( Kontar. Hannah, and MacKinnon . 2008 ; Kontar et aZ.I . [2010h . Kontar et al. an- 
alyzed a GOES M6 class flare observed by RHESSI on the east solar limb. They 
found observed source heights from 1 to 0.5 Mm for source from 20 to 200 keV 
consistent with thick-target emission. 

Here, we consider the RHESSI microflare height distribution. We use the 
RHESSI microflare list which is de scribed and was prev iously analyzed in two 
papers bv lChriste et all (|2008l) and lHannah et all (|2008l) . 
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Figure 1. Three examples of flare position finding. Positions are found from back-projection 
images of each event integrated over 4 minutes around the peak of the flare (in the 6—12 keV 
band). Images were created using detectors 5, 6, 8, and 9 individually (detector 7 has little 
sensitivity at low energies). For each image, a paraboloid was fit to the maximum in the image. 
Here, individual detector positions are shown overlaid on a representative back-projection 
image with 60%, 70%, 80%, 90% contours. (Left) An unresolved flare where positions from 
each detector agree. (Center) An unresolved flare where the flare position from detector 9 
disagrees with the others due to the flare's proximity to the spin axis. (Right) A resolved 
flare where positions from each detector disagree with each other since each detector resolves 
different sizes of the source. 



2. Observations 



RHESSI consists of 9 ultra- pure germanium detectors ( Smith et ai , 20021). each 
behin d a rotation modulation collimator (RMC) which enable imaging ( Hurford et q/J .[ 
20021 ). Each RMC is made up of a pair of identical grids of slits and slats with 
different pitches which modulate incoming X-rays as the spacecraft rotates. Each 
RMC measures a different spatial frequency down to 2.3 arcscconds. 

RHESSI has been observing the Sun since its launch in February 2002. A 
central product of the RHESSI collaboration is the so-called RHESS I flare list. 
A separate RHESSI microflare list also exists ( Christe et all . 120081) . This list 
contains the start, peak, end times, and background times along with positions of 
all microflares observed by RHESSI from March 2002 to March 2007. It contains 
25,006 events of which 94.3% have associated positions as found in the 6-12 keV 
energy range. Small flare events with simple signatures of nonthermal electrons 
and thermal heating are well suited to study elementary flare processes compared 
to larger events with complex morphologies. 

Here, positions were redetermined for each microflare. Positions were found 
through back projection images of each event integrated over 4 minutes around 
the peak time in the energy ranges, 4-10, 10-15, and 15-30 keV. The median 
nonthermal low energ y cutoff for microflares was previously found to be 12 keV 
([Hannah et all . 120081 ). therefore the lowest energy band is dominated by thermal 
emission while the highest energy range is likely nonthermal. Images were created 
using detectors 5, 6, 8, and 9 individually (detector 7 has little sensitivity at 
low energies and was excluded). The angular resolutio ns are 20.4, 35.3, 1 05.8, 



183.2 arcseconds, for detectors, 5, 6, 8, 9, respectively ([Smith et ai , l2002t) . For 



each image, a paraboloid was fit to the maximum in the image returning an 
interpolated flare position with errors (see Figure Q] for examples). 

The median flare loop length and w idth for micro f lares have been found to 
be 32 and 11 arcseconds, respectively ([Hannah et al ., 2008). This implies that 
detectors 6 and above may not resolve flare loops and should provide centroid 
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Figure 2. The locations of RHESSI microflares from the RHESSI microflarc list 
(Christe et al., 2008) rc-cvaluatcd using back projections images from detectors 5, 6, and 8 
and averaging the positions from each. Those events with a standard deviation in position 
greater than 5 arcseconds were excluded. Flare positions are concentrated in the active region 
bands and near the limb. A lack of events can clearly be seen in the 4-10 keV positions 
around [250,— 100] arcsec, the average location of the spin axis, where imaging is difficult 
dChriste et all 120081 '). 



positions that agree with each other. This fact provides a necessary cross check 
since it is not straightforward to algorithmically determine whether centroid lo- 
cations which disagree between detectors are due to a resolved source or a noisy 
image. In order to include only those flares with confirmed positions, results from 
detectors 5, 6, 8, and 9 were compared. Those events with a standard deviation 
in position greater than 5 arcseconds (half of the median flare loop width) were 
excluded. This reduced the number of events from 25,006 to 8,503 for the 4- 
10 keV energy range. Detector 9 positions were found to differ significantly from 
other detectors. This is likely due to the fact that detector 9 has the coarsest 
resolution and is therefore most often affected by the spin axis (where imaging 
i s difficult). The ave rage location of the spin axis is around [250,-100] arcsec 



(jChriste et all 120081 ). therefore flares on the west limb are 4 times the resolution 



of detector 9 away from the spin axis and should be largely unaffected. The final 
dataset now using detectors 5, 6, and 8 includes 14,750 events for 4-10 keV, 
5,829 for 10-15 keV, and 1,452 for 15-30 keV. The flare positions can be seen in 
Figure [51 Errors in the fitted centroids are small for all detectors. The average 
error in flare positions for detector 5 is 0.68 arcsec. These errors are smaller than 
the difference in positions between detectors, therefore the standard deviation 
of the set of positions is used as the error instead. 



3. Discussion 



Observations of solar source positions are a convolution of the radial source 
height and the spherical shape of the Sun. It is therefore difficult to determine 
the absolute height above the photosphere for individual events. In spherical co- 
ordinates, the relevant variables are the azimuthal angle, 9 (or 77, the longitude), 
the polar angle, 4> (<fi = ir/2 — 8 where 6 is latitude), and the radius, r, which is 
related to the source height, h, by r = Rq + h where Rq is the solar radius. The 
position of an event is fully described by (r, 77, </>). At the Earth, the only ob- 
servable quantity is the two dimensional projection (X = (Rq + h) sin(?7) cos(r5), 
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Figure 3. The radial distribution of observed flare positions near the limb for each energy 
range (4-10 kcV, 10-15 keV, 15-30 keV). The positions were corrected for the apparent size 
of the Sun and the solar B angle. The dashed line represents the optical solar limb. Positions 
are found to be concentrated at the limb with a peak at a heliocentric radius of approximately 
1.OO5-R0 and a steeply falling distribution above the limb. The peak in the distribution above 
the limb is an artifact due to partially-occulted sources. 



Y = (Rq + h) sin(<5)) of this vector. The height of a source is easiest to observe 
at the limb where projection effects arc small. In order to isolate the limb, we 
first define the radial distance of the source from the center of the solar disk as 
the projected source radius, p = y ' X 2 + Y 2 , which is related to the coordinate 
variables through the following equation, 



p = {R Q + h) ^/sin^) 2 cos(<5) 2 + sin(,5) 2 . (1) 

The measured projected source radius distribution for each energy band can be 
seen in Figure [3J Most events are observed to occur above the limb. 

Generally none of the variables on the right-hand side of equationQ]are known. 
Yet if the distributions for each variable (g{T]),g(S),g(h)) is known then the 
distribution of the projected source radius, g(p), can be determined. Assuming 
all events are on the equator (g(6) = 0), have identically zero height (g(h) = 0), 
and are uniformly distributed in longitude, the distribution function for p is 
given by 

9(P) = ; A (2) 

where A is a normalization constant. Generally, deriving an analytical solution 
for g(p) is impossible, therefore a forward fit model was developed. The model 
takes as input a source height distribution, g(h), a longitude and latitude dis- 
tribution and generates random vectors (h, 77, S). The longitude distribution 
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is assumed to be uniform^ while random la titudes are generate d based on the 
observed microflare latitude distribution (see lChriste et aZ.l . l2008l ). These random 
vectors are then projected onto the two dimensional visible solar disk as (X, 
Y) excluding those that are occulted. The source projected radius, p, is then 
calculated for each event and the distribution, g(p), can be compared with the 
measured distribution to determine the height distribution, g(h). 

In this study, two height model distributions are considered; an exponential 
distribution defined by g{h) = A exp(^) and a power-law distribution, g{h) = 
Ah~ a . These t wo models were chosen for their mathematical simplicity following 
the example of lCatalano and van Alien! (|l973h and are not necessarily physically 
motivated. 

The discussion so far has assumed point sources. Finite source size should 
have no effect on the projected radial distribution for sources well above the 
photosphere (p > a size)- Those sources whose maxima are behind the limb yet 
are within ~cr S j Ze will show emission truncated by the solar limb. RHESSI imag- 
ing of a partially occulted source will be "blurred" and show a source maximum 
above the limb whose height is dependent on the detector resolution. This effect 
places a minimum threshold below which the height of occulted sources cannot 
be accurately determined by RHESSI and is the reason for the maximum in the 
observed radial distribution. For detector 5, 6, and, 8 the height of the peak in 
the radial distribution was found to be 5.0, 5.3, and 6.7 arcsec above the limb, 
respectively. Height determination for sources as close to the limb as these values 
should be investigated carefully in order to determine whether they are occulted. 
In the following analysis, both height models include an unphysical minimum 
cutoff height, ho, to model this effect. Using only unresolved sources, as is done 
here, is unlikely to affect the shape of the height distribution though there may 
be an affect on the minimum source height value if the source height is related 
to the source size. 



4. Analysis 

Using the method described above, the height model distributions were fit to 
the observed distribution through multiple Monte Carlo runs. Best fit parame- 
ters were determined through a grid search for the minimum in the x 2 space. 
The number of simulated events required for the Monte Carlo simulation were 
chosen in order to make the variance in the minimum smaller than the errors 
as determined by the contours of the x 2_s P ace while maintaining an accept- 
able computation time. For the exponential height model, in the 4-10 keV 
energy range, the scale height, A, was found to be 6.1 ± 0.3 Mm and the min- 
imum height, ho, is 5.7 ± 0.3 with a reduced \ 2 °f 1-06. The fit results can 
be seen in Figure 2] The authors would like to point out the presence of a 



It is incorrect to select longitude from a uniform distribution in order to have points dis- 
tributed uniformly over a sphere since the spherical area element depends on latitude. This 
effect is small (up to a few percent) since most flares occur near the solar equator and 
simulations show that it does not affect the results presented here. 
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Figure 4. The radial distribution of observed flare positions near the limb in the 4—10 keV 
(black) energy bands. The dashed line represents the optical solar limb. (Left) The best fit 
assuming an exponential height distribution is overlaid in blue. The scale and minimum height 
were found to be 6.1±0.3 Mm and 5.7±0.3 Mm, respectively. The exponential fit captures the 
overall behavior of the data though underestimates the number of events at large heliocentric 
radii. (Right) The best fit assuming a power-law height distribution. The power-law index and 
minimum height were found to be 3.1 ±0.1 Mm and 6.9 ±0.3 Mm, respectively. The power-law 
height distribution fit is worse than the exponential height distribution as measured by the 
reduced x 2 but better captures the events at large heliocentric radii. 



step-like structure at w 0.96i?©. This feature is statistically significant and is 
most prominent in the eastern hemisphere. Its cause is unknown to the au- 
thors but it does n ot affect the results presented h ere. The fit results disagree 
with those found bv lCatalano and van Allen! ( 19731) who found a scale he i ght o f 
11 ±3 Mm in the 1-6 keV range thoueh lAschwanden. Brown, and Kontarl (|2002n 
concluded that most past height m easurements have large systema tic systematic 
errors. Ha flare positions, which ICatalano and van Allenl (j 19731) depends on, 
have an ass ociated accuracy of 1 he liographic degree (or 12 Mm). Analyzing 
the study by Matsush ita et al. ( 1992), w hich similarly depend on Ha positions, 
Aschwanden. Brown, and Kontarl"( 20021) requi r ed a c orrection factor of —6.8 Mm 



which, if applied to ICatalano and van Allenl (|l973l) . would lead to agreement 
with the results found here. 

Fits were also performed for the 10-15 keV and 15-30 keV energy range. 
Though those results have poorer statistics, they agree with the 4-10 keV fit to 
within la. Thes e results are consi stent with past observation s of thermal loops 
such as those of Ohki et al. ( 1983 ) and Takakura et al. ( 19861 ). 

The power-law height distribution was also fit to the data for each energy 
channel. The power-law index, a, was found to be —3.1 ± 0.1 with a minimum 
height, ho, of 6.9 ± 0.3 Mm, similar to the minimum height found in the cxpo- 
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Figure 5. The displacement versus projected radial position for each energy range. A fit to 
the data yields the average height difference between energy ranges (white line) . The errors bars 
represent the standard deviation of the positions from different detectors. The average height 
differences are small, A/i4_io,io— 15 keV = f 50±25 km, Ahio_i5 15—30 keV = 270±51 km, and 
A/i4_io i5_3o keV = 320 ± 55 km where lower energies are at larger heights. The results arc 
found to be self- consistent since A7i 4 _i 16 _ 30 kc v = ^4-10,10-15 keV + Afeio-15, 15-30 keV- 



ncntial height model, and a reduced \ 2 of 1.31. Again, the results at the higher 
energies were found to agree. Both the exponential and power-law height models 
fit the data well with a small reduced \ 2 - The exponential fit was found to fit 
the data close to the limb better than the power-law model. This fact explains 
the lower value of x 2 for the exponential fit. On the other hand, the power-law 
fit better captures the tail of events at high heliocentric radius. 

To better compare the sour c e heig hts between different energies, a method 
develo ped bv iMatsushita et ~al. (1992) and more recently applied bv lSato et al 



( 20061 ) is appropriate here. Assuming that multiple pairs of points lie on either 
one of two concentric spheres with different radii, then the relationship between 
the projected radial positions and the length of the displacement vectors between 
each pair of points is related to the difference in radius (or height, h) between 
the two spheres. This method, by design, ignores those events observed above 
the limb and is therefore not sensitive to the problem of partially occulted 
sources. Figure \5\ displays this relationship for each energy range. The errors bars 
represent the standard deviation of the positions from different detectors. The 
fitted height differences were found to be small; A/i4_io,io-i5 koV = 150±25 km, 
A/iio-15, 15-30 koV = 270 ± 51 km, and Aft, 4 _io, 15-30 kcv = 320 ± 55 km where 
lower energies are at larger heights. These results are found to be self-consistent 
since it must be that A/i 4 _i , 15-30 kcv = A/i 4 _ 10,10- 15 kcv + A/iio-15, 15-30 kcv- 
The reduced x 2 f° r ai l of the fits are less than one which suggests that the errors, 
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Figure 6. The height distribution models compared. The fitted exponential (black) and 
power-law (red) models (a = 3.1 ± 0.1) are shown. The two models diverge significantly 
above O.O5i?0. The PFSS loop height distribution is also shown along with a power-law fit 
(otpFSS = —1-4 ± 0.12). The PFSS loop distribution has a much harder distribution than the 
measured distribution which may be due to the energization process of the corona. 



which here represent the standard deviation in the set of detector positions, are 
too large. This may indicate that source structure is being resolved. Comparing 
to previous uses of this met hod, it is found that the heig ht differences are small 
compared to results by both Matsushita et al. ( 19921 ) and Sato et al. (|2006l) . The 
absolute heights in each energy band suggests that images in all energy bands 
are dominated by thermal sources. 

In order to better interpret the fitted height distribution, we used the potential-| 
field source-surface (PFSS) to model (non-flaring) magnetic loops. This model 
calculates the magnetic field in the corona using photospheric magnetic field 
observations (SOHO/MDT) assuming that there are no currents (V x B = 
0). This model was initially deve loped bv lAltschuler and Newkirk (1969) and 
Schatten. Wilcox, and Nesa ( 19691 ) , and later refined by Hoeksema ( 1984 ) and 
Wang and Sheelevl (jl992l )~The PFSS model, through its assumption of zero cur- 
rent, models the lowest energy state of the (non-flaring) corona. It is, therefore, 
not an accurate model of the actual corona below ~1.6 solar radii yet it provides 
a useful theoretical "base" corona to compare with. 



Using the standard PFSS softwar^l developed by M. DeRosa (see Schriiver and De RosaL 



20031 ) and distributed with SolarSoft, the PFSS height model was calculated from 
0.01 to 2.Oi?0. In order to get an accurate representative corona, the PFSS model 
was calculated every thirty days from March 2005 to March 2007, the time range 
of the RHESSI microflare list. Magnetic field lines were traced for each date and 
the maximum height for every (closed) loop was determined. From this the PFSS 
loop height distribution was then determined and each was fit with a power-law. 



2 see http://www.lmsal.com /~derosa / pfsspack 
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The distribution shown here is the average for every date. The power-law index, 
ctpFSS: °f the average distribution was found to be —1.4 ± 0.12. The error in 
the power-law index represents the standard deviation of the set of fit values 
for individual dates. The distribution of PFSS loop-top heights can be seen in 
Figure [5] compared with the power-law and exponential fits to the microflare 
height distribution. The PFSS height distribution was found to be much less 
steep than the power-law fit to the data with a difference in power-law index of 
1.8 ± 0.16 which implies that the ensemble of potential loops contain more large 
loops than the observed distribution. A soft loop height distribution may be a 
signature of the energization process of the corona. 

5. Conclusions 

Taking advantage of the large dataset provided by the RHESSI microflare list, 
the height distribution of HXR microflares was determined. The distribution 
of HXR (4-10 keV) source heights was found to be well fit by an exponential 
distribution with a scale height of 6.1 ±0.3 Mm. A power-law height distribution 
with a negative power-law index of 3.1 ±0.1 was also found to be consistent with 
the data though a worse fit. The minimum observable height due to partially 
occulted sources was found to be 5.7 ± 0.3 Mm. The value of the heights found 
here suggests that they are thermal loop sources. Thermal loops are frequently 
assumed to be semicircular loops. The center of mass of a semicircular loop 
is located 36% below the loop-top therefore the minimum observable loop-top 
height may be as low as 7.7 Mm. In practice, RHESSI images of loops do not 
usually show the base of the loop, therefore it is likely that this is an overestimate. 

The height distributions were compared to the loop height distribution as 
determined by a potential field model (PFSS). It was found that the measured 
flare height distribution was much steeper than the PFSS magnetic loop heights 
(difference in the power-law index of ~1.8) meaning the ensemble of potential 
loops contain many more large loops compared to the height distribution found 
here. Observations of rising post-flare loops in EUV suggest that the flare loop 
distribution may relax to the PFSS distribution as energy is dissipated. In order 
for the observed distribution to transition to the potential distribution small 
flare loops must rise faster than large flare loops. It must also be true that an 
energized corona suppresses large loops. 

The height distribution for different energies (4-10 keV, 10-15 keV, 15- 
30 keV) were found to agree within statistical uncertainties. More sensitive analy- 
sis showed that the difference in the average height were small, Aft,4_io,io-i5 keV =| 
150±25 km, A/iio-15. 15-30 keV = 270±51 km, where lower energies are at greater 
heights. This suggests that emission in each energy range should be interpreted 
as from thermal flare loops. 

The minimum observable height found, = 6.3 ± 0.3 Mm, introduced by 
partially-occulted flares, suggests that investigations of flares at the limb should 
be mindful of this limitation whereby limb sources may appear at higher altitudes 
if they are partially behind the limb. Such events should be detectable by the fact 
that images with different detectors will show positions which move away from 
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the limb as a function of resolution. It should also be possible to detect partially- 
occulted flares through the visibilities in the highest resolution detectors which 
should have unusually large amplitudes for directions parallel to the limb. 

Acknowledgements This work was partially supported under NASA grant NAS5-98033 
and NNM05ZA12H. The authors would like to thank the referee for input which has greatly 
improved this paper. 

References 

Altschulcr, M.D., Newkirk, G.: 1969, Solar Phys. 9, 131. 

Aschwandcn, M.J., Brown, J.C., Kontar, E.P.: 2002, Solar Phys. 210, 383. 

Aschwanden, M.J., Schwartz, R.A., Alt, D.M.: 1995, Astrophys. J. 447, 923. 

Catalano, CP., van Allen, J. A.: 1973, Astrophys. J. 185, 335. 

Christc, S., Hannah, I.G., Kruckcr, S., McTiernan, J., Lin, R.P.: 2008, Astro- 
phys. J. 677, 1385. 

Hannah, I.G., Christe, S., Krucker, S., Hurford, G.J., Hudson, H.S., Lin, R.P.: 
2008, Astrophys. J. 677, 704. 

Hocksema, J.T.: 1984, Structure and evolution of the large scale solar and 
heliosphcric magnetic fields. PhD thesis, Stanford Univ., CA.. 

Hurford, G.J., Schmahl, E.J., Schwartz, R.A., Conway, A. J., Aschwanden, M.J., 
Csillaghy, A., Dennis, B.R., Johns-Krull, C, Kruckcr, S., Lin, R.P., Mc- 
Tiernan, J., Metcalf, T.R., Sato, J., Smith, D.M.: 2002, Solar Phys. 210, 
61. 

Kane, S.R.: 1983, Solar Phys. 86, 355. 

Kane, S.R., Anderson, K.A., Evans, W.D., Klebesadel, R.W., Laros, J.: 1979, 
Astrophys. J. Lett. 233, 151. 

Kontar, E.P., Hannah, I.G., MacKinnon, A.L.: 2008, Astron. Astrophys. 489, 
57. 

Kontar, E.P, Hannah, I.G., Jeffrey, N.L.S., Battaglia, M.: 2010, Astrophys. J. 
717, 250. 

Matsushita, K., Masuda, S., Kosugi, T., Inda, M., Yaji, K.: 1992, Pub. Astron. 
Soc. Japan 44, 89. 

Ohki, K., Takakura, T., Tsuneta, S., Nitta, N.: 1983, Solar Phys. 86, 301. 

Sato, J., Matsumoto, Y., Yoshimura, K., Kubo, S., Kotoku, J., Masuda, S., 
Sawa, M., Suga, K., Yoshimori, M., Kosugi, T., Watanabe, T.: 2006, Solar 
Phys. 236, 351. 



SOLA: paper.tex; 19 January 2013; 17:52; p. 11 



Christe et al. 



Schatten, K.H., Wilcox, J.M., Ness, N.F.: 1969, Solar Phys. 6, 442. 

Schrijver, C.J., Dc Rosa, M.L.: 2003, Solar Phys. 212, 165. 

Smith, D.M., Lin, R.P., Turin, P., Curtis, D.W., Primbsch, J.H., Campbell, R.D., 
Abiad, R., Schroeder, P., Cork, CP., Hull, E.L., Landis, D.A., Madden, N.W., 
Malone, D., Pelil, R.H., Raudorf, T., Sangsingkeow, P., Boyle, R., Banks, I.S., 
Shirey, K., Schwartz, R.: 2002, Solar Phys. 210, 33. 

Takakura, T., Tanaka, K., Nitta, N., Kai, K., Ohki, K.: 1986, Solar Phys. 107, 
109. 

Tsuneta, S., Takakura, T., Nitta, N., Makishima, K., Murakami, T., Oda, M., 
Ogawara, Y., Kondo, I., Ohki, K., Tanaka, K.: 1984, Astrophys. J. 280, 887. 

Wang, Y., Shcclcy, N.R. Jr.: 1992, Astrophys. J. 392, 310. 



SOLA: paper.tex; 19 January 2013; 17:52; p. 12 



